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We apply the effective field theory approach to quasi-single field inflation, which contains an addi- 
tional scalar field with Hubble scale mass other than inflaton. Based on the time-dependent spatial 
diffeomorphism, which is not broken by the time-dependent background evolution, the most generic 
action of quasi-single field inflation is constructed up to third order fluctuations. Using the obtained 
action, the effects of the additional massive scalar field on the primordial curvature perturbations are 
discussed. In particular, we calculate the power spectrum and discuss the momentum-dependence 
of three point functions in the squeezed limit for general settings of quasi-single field inflation. Our 
framework can be also applied to inflation models with heavy particles. We make a qualitative 
£Nj ■ discussion on the effects of heavy particles during inflation and that of sharp turning trajectory in 

our framework. 
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Inflation gives the most natural solution to the horizon and the flatness problems of the big-bang theory as well as 
generates the primordial perturbations [J , whose properties well coincide with the recent observations of cosmic 
microwave background anisotropics like the Wilkinson Microwave Anisotropy Probe Models of inflation can be 
classified into two categories with respect to relevant degrees of freedom during inflation: single-field models and 
Q_[ multiple field models. Recently, the most general single field inflation model with the second order equations of 
motion [|J has been invented in the context of Horndeski 0, [f| and Galileon theories @, Q ■ Then, the bispectra of 
i ' primordial curvature 0, [loj and tensor perturbations [Tlj are obtained as well as their powerspectra [|J . 

\ Effective scalar fields are ubiquitous in the extensions of the standard model of particle physics such as supergravity 



> 



and superstring. Then, it is well motivated to consider multiple field models of inflation. Such multiple field models 
are roughly divided into three classes: (i) only one field is light, while the other fields arc very heavy compared to the 
(N| . Hubble scale during inflation. Generically, this class virtually falls into the single field category fl2ll . However, it is 
' recently discussed that heavy modes can affect the dynamics of light mode in some particular cases (ii) there 

are multiple light fields, in which isocurvaturc perturbations arc generated in addition to curvature perturbations, 
(iii) only one field is light, while the masses of other fields are comparable to the Hubble scale during inflation. This 
class is called quasi-single field inflation model [2l|, [22[ . 

(N- 

In supergravity, inflation necessarily involves sup ersymmetry (SUSY) breaking, whose effects are transmitted into 
other scalar fields as Hubble induced masses [23M25I ] . 1 Therefore, quasi-single field inflation is naturally realized 
in supergravity and it is well motivated by the model building based on supergravity or inspired by superstring. 
■ Furthermore, it is known that massive isocurvature modes which couple to the inflaton and have Hubble scale masses 
$H \ can give significant impacts on primordial curvature perturbations. In the original paper plj by Chen and Wang, it 
was shown that, for example, scalar three point functions take the intermediate shapes between local and equilateral 
types. Based on these backgrounds, we would like to discuss quasi-single field inflation model in general settings. 



Recently, effective field theory approach to inflation has been invented in |27h29I |. which is based on the symmetry 
breaking during inflation: Time diffeomorphism is broken by the time-dependent background evolution during infla- 
tion. Then, based on the unbroken time-dependent spatial diffeomorphism, the effective action for inflation can be 
constructed systematically in unitary gauge, where inflaton is eaten by graviton and there are no perturbations of 
inflaton. By use of the Stiickclbcrg trick, the curvature perturbation can be associated with the Goldstone boson 7t, 
which non-linearly realizes the time diffeomorphism. The key observation is that the Goldstone n could decouple 
from the metric fluctuations in some parameter region which we call the decoupling regime. In the decoupling regime, 
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1 The methods to keep inflaton flat against such SUSY breaking effects are reviewed in Refs. [2^1 . 
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the dynamics of Goldstone 7r is described by a simplified action, which does not contain metric perturbations. As a 
consequence, calculations of scalar perturbations are also simplified and seeds of non-Gaussianities become clear. 

In this paper, we apply this effective field theory approach to quasi-single field inflation. First, in unitary gauge, 
we write down the most general action invariant under the time-dependent spatial diffcomorphism and constructed 
from graviton and the massive isocurvature mode. The obtained action is expanded systematically in fluctuations 
and derivatives around the FRW background. By the Stuckclbcrg trick, wc introduce the the action for Goldstone 
boson and carefully discuss its decoupling regime. Using the action in the decoupling regime, the power spectrum 
is calculated in the general setting of quasi-single field inflation. The momentum dependence of scalar three-point 
function is also discussed in the general setting. Our framework can be also applied to inflation models with heavy 
particles. As an application, we make a qualitative discussion on the effects of heavy particles during inflation and 
that of sharp turning trajectory. 

The organization of this paper is as follows. In the next section, we briefly review the effective field theory approach 
and quasi-single field inflation. In Scc. lIIIl the most general action for quasi-single field inflation is constructed via 
effective field theory approach. The decoupling regime of the obtained action is also discussed. In Sec. HVl the power 
spectrum is calculated first in the general setting of quasi-single field inflation with constant mixing couplings. Then, 
the effects of sharp turning trajectory on the power spectrum is qualitatively discussed. In Scc.|Vl the momentum 
dependence of scalar three-point functions arc discussed. Final section is devoted to summary and discussions. 
Technical details of the calculation of the power spectrum are summarized in Appendices. 



II. EFFECTIVE FIELD THEORY APPROACH AND QUASI-SINGLE FIELD INFLATION 

In this section we briefly review the effective field theory approach to inflation developed in [27| and the quasi-single 
field inflation model proposed in (2l| . 



A. Effective field theory approach to inflation 



Inflation is an accelerated cosmic expansion with an approximately constant Hubble parameter: 

ds 2 = -dt 2 + a 2 {t)dx l dx l with H(t) = - , e = --^-<l. (1) 

a H 2 

It is characterized by the spontaneous breaking of the timc-diffcomorphism: 

=M*), (2) 

where <p(t, x) is some scalar operator. Here we chose the frame in which the vacuum expectation value of <j)(t, x) is 
spatially uniform. Assuming the degrees of freedom relevant to the cosmological perturbation and invariancc under 
the time-dependent spatial diffcomorphism, x l — > x' 1 = x l + £*(f, x 3 '), which is not broken by the condensation <fio(t), 
we can construct the effective action for inflation. 

In the simplest case, relevant degrees of freedom are three physical modes of graviton: two transverse modes and 



one longitudinal mode related to the inflaton. As discussed in |27| , any action of graviton invariant under the time- 
dependent spatial diffcomorphism can be written in terms of the Ricmann tensor R^pa, the time-like component 
of the metric g 00 , the extrinsic curvature on constant-^ surfaces, the covariant derivative V M , and the time 
coordinate t: 

S = J d 4 x^gF(R^ pa ,g 00 ,K^,Vp,t), (3) 

where all the free indices inside the function F must be upper O's. Note that g 00 should be treated as a scalar when 
considering its covariant derivative, and we can use d^g 00 for example. The explicit form of the extrinsic curvature 
is 



<%0„ g ° o + ff0 Mg ° Q ffffg 0g a g g 00 g ?(d,g p „+d„g p »-d p g^) 
P^iv -neatly 2(-g 00 ) 3 / 2 2(- 5 00 ) 5 /2 2{-g 00 ) 1 / 2 
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where n p = — —=!=== is a unit vector perpendicular to constant t surfaces and h pv = g^ v + n^,n v is the induced 

spatial metric on constant t surfaces. In [27], it was shown that the action ([3]) can be expanded around a given FRW 
background as 



5 



d x-s/—g 



l-M&R + M&Hg 00 - M 2 X (3H 2 + H) + F® '(Sg 00 \5K pu ,8R pvp<J - S^g^g^ V ^t) 



(5) 



where the function starts with quadratic terms of the arguments Sg 00 , SK pu , and 5R pypa and all the free indices 
must be upper O's. The arguments Sg 00 , 5K pl ,, and 8R pvpa are defined by 2 



Sg 00 = g a 



1 



5K pv = K pu - Hh pu , 
SR^pa = R^pa - 2H 2 h p[p h a ] v + (H + H 2 ){h pp 5 v 5 + (3 permutations)) 



(0) 
(7) 
(8) 



They are covariant under time-dependent spatial-diffeomorphism and vanish on the FRW background. Notice that 
the action of single field inflation in the uniform inflaton gauge can be reproduced by gauge-fixing the time-dependent 
spatial diffcomorphism as 

g ij (x)=a 2 (t)e 2 ^(e^) ij with 7ii = £> i7 y = , (9) 

where Ci x ) is the scalar perturbation. 

For the calculation of correlation functions of the scalar perturbation f , it is convenient to introduce the action for 
the Goldstone boson ir by the Stiickelberg method. We perform the following time-diffcomorphism on the action (JSJ) 
in the unitary gauge: 



t ->• i , x 1 -> x % with i + 7r(t, x) = t , x l = x % . 
In general, the transformation (flTj|) is realized by the following replacement: 



(10) 



V 



c? x\/—g , 



i? 



fj,vpcr 



R 



where we dropped the tilde for simplicity and g transforms, for example, as 
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g 00 + 2g°^d p n + g^d p nd u n. 



(11) 



(12) 



The transformation rules of K pv and h pv also follow from (|TT|) straightforwardly. These procedures lead to the 
following action for the Goldstone boson ix: 



S 



d X\J—g 



-M 2 X R + M 2 x H{t + tt) (g aa + 2g 0fl d p ir + g^d^d^) - M 2 Y ('iH 2 {t + tt) + H(t + tt) 



(13) 



where the dots stand for the terms corresponding to F^ 2 \ The obtained action enjoys the timc-diffeomorphism by 
assigning to tt the non-linear transformation rule 3 



tt(.t) -> tt(.t) = tt(x) - with i ^ i = i + £°(x) , x i 



(14) 



and the action in the unitary gauge can be reproduced by gauge-fixing the time-diffeomorphism as ir(x) = 0. 



2 Here and in what follows, we concentrate on the spatially flat FRW background. 

3 Here it should be noted that while the action is invariant under the time-diffeomorphism, it does not have the shift symmetry 
tt — > tt + constant because of the time-dependent free parameters such as H(t + tt) or H(t + tt). Expanding the parameters in it, we 
find, for example, that tt has a mass term ~ Mpjif 2 7r 2 , which is sub-leading in the slow-roll approximation. 
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It is important to recognize that, in the action (|13[) . terms with graviton fluctuations have less derivatives than 
those without graviton. Because of this property, it is expected that the mixing of the Goldstone boson and graviton 
becomes irrelevant to the dynamics of the Goldstone boson at a sufficiently high energy scale. For example, let us 
consider the following simplest case: 



S 



d X\J — g 



-M^R + M^H(t + tt) (g 00 + 2g "d ll w + g^d^d^) 



Ml, 



(3H 2 (t 



H(t 



(15) 



In the canonical normalization (tt c ' 
written as 



Mpii-H) 1 / 2 **, 8g% v ~ Mpidg^), the mixing term M*H 5g°>* 9„7T can be 



(16) 



and it can be neglected in the energy scale E 3> e 1 ^ 2 !!. In other words, when the slow-roll parameter e is small, 
the mixing becomes irrelevant inside the horizon. The dynamics of tt inside the horizon are then determined by the 
following action in the decoupling limit: 

(^tt) 2 * 



M^H{t + TT)[ -1-27T-7H 



- M^(3H z (t + tt) + H(t + tt)) 



(17) 



where the metric reduces to the FRW background. More generally, the mixing of the Goldstone boson and graviton 
becomes irrelevant inside the horizon when the free parameters of the action are in some regime (decoupling regime) 
as in the slow-roll regime for the above simplest case. In the decoupling regime, the calculation of scalar correlation 
functions also becomes tractable. Taking the spatially fiat gauge 4 



g i:j (x) = a 2 {t){e l{ - x) ) tj with lH = d^tj = , 



(18) 



the scalar perturbation £(x) is given by £(x) = —H{t)Tr{x) at the linear order, and the calculation of correlation 
functions of £ reduces to those of tt, which can be obtained using the simplified action in the decoupling limit. This 
kind of simplification in the decoupling regime is one of the advantages to use the effective field theory approach. In 
the next section we extend this approach to quasi-single field inflation. 



B. Quasi-single field inflation 

The original model [2l[ of quasi-single field inflation is described by the following matter action: 

^matter " 



d 4 xV=g[ - \{R + X?9^d^9d v e - \y»d^ X dvX - V„{9) - V(x) 



(19) 



where 6 and x are the tangential and radial directions of a circle with radius R and the potential V sr (6) along the 
tangential direction is of slow-roll type. The homogeneous backgrounds and their equations of motion are given by 

9 = 6 (t), x = Xo (constant) , 

3M 2 ,H 2 = ±R 2 2 + V( Xo ) + V ST (Oo) , -2M 2 ,H = R 2 9 2 , V'( X o) = RB% , R% + 3R 2 H9 Q + Vj T (9 ) = , (20) 

where R = R + xo- Expanding the action around the homogeneous background, it yields the following second order 
action of the fluctuations 59 = 9 — 9 and a = x ~ Xo : 



,(2) 



l -R 2 g^d^59d u 59 - l -g^d^od v o - R9 ad°S9 - ^(V"( X o) - 2 )a 2 



(21) 



The mixing coupling aS9 converts the er 3 coupling, for example, into three point functions of 69, and hence this model 
can potentially give a large non-Gaussianities. Furthermore, it is known that the squeezed limit of scalar three point 
functions is sensitive to the mass of a: 



lim (CkiCkaCks) a K 3/2 " k i 6 > 

K3/k 1 =k 3 /fc 2 =K— >0 



(22) 



4 Strictly speaking, the name of this gauge may be inadequate because there are still tensor fluctuations and hence the spatial hypersurface 
is not exactly flat. 
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where v = y — — and m| = V"(xo) — R@o < 1^- As this snTL pl e model implies, massive scalar fields of Hubble 

scale mass can cause a non-trivial behavior of non-Gaussianities. In the following sections we discuss more general 
setting for quasi-single field inflation using the effective field theory approach. 



III. MOST GENERIC ACTION OF QUASI-SINGLE FIELD INFLATION 

In this section, we construct the most generic action of quasi-single field inflation using the effective field theory 
approach. After constructing the action in the unitary gauge first, we derive the action for the Goldstone boson tt 
and discuss its decoupling regime. Relations between our approach and models in the literatures are also discussed. 



A. Action in the unitary gauge 



In the unitary gauge, the relevant degrees of freedom in quasi-single field inflation are three physical modes of 
graviton and an additional scalar field a. The typical mass of a is supposed to be of the order of the Hubble 
scale during the inflationary era. In this subsection, we construct the most generic action invariant under the time- 
dependent spatial diffcomorphism from graviton and the scalar field a up to the third order fluctuations. Here it 
should be noticed that the action constructed in this section can be applied to any two field models because no 
conditions on a arc imposed. 5 

Extending the procedures in [27j to our case, the most general action invariant under the time-dependent spatial 
diffcomorphism is given by 



d 4 xy^g F(R pvpa ,g 



and it is expanded around the given FRW background as 



do 



S = / d 4 xJ=l 



\m^R + M^Hg 00 - M^(3H 2 + H) + F^(5g QQ ,a, SK^SR^X 



(23) 



(24) 



where all the free indices inside the functions F and F^ must be again upper O's and F^ starts with quadratic 
terms of the arguments Sg 00 7 a, SK^, and SR^ pa . 6 Then, let us write down possible terms in the action up to the 
third order fluctuations. Schematically, we write the action in the following way: 



S — Si 



grav 



S(j + Sri 



(26) 



where the first term S glav in the right hand side denotes terms constructed from 5g , 5K^ U , and 5R^ vpa , the second 
term S a denotes those only from a, and the last term S'mix denotes those mixing the graviton fluctuations and a. As 
discussed in [27j, the first term 5 grav can be expanded as 



d 4 Xy 



^M^R + M^H(t)g 00 - Mfi (^3H 2 (t) + H{t) 



MM 
2! 



00\3 



(V U ) 



(27) 



5 For example, we do not require the shift symmetry of cr. cr — > cr + constant, which is assumed in multi-field inflation [28l . The mass 
of cr is not necessarily of order Hubble scale so that the action constructed in this section can be applied not only for quasi-single field 
inflation but also inflation models with an additional heavy scalar. 

6 In general, sums of terms linear in the fluctuations can be practically second order. For example, let us consider the term 
J d 4 Xy/—g [/i(i)cr + f2(t)d°a\. Although this kind of action seems to be first order apparently, it turns out to be second order af- 
ter taking into account the equation of motion for cr: fi(t) + f2(t) + 3Hf2(t) = 0. Then, the function F( 2 ) seems to contain such a 
combination of linear order terms. However, using the relation 

j d 4 Xv^/(t)V°(. •■) = -/ «J 4 aV=ff/(t)V Z i Bff n' i S#.C ■ •) = / <* 4 *V=s( - g m f(t) - \f{t)d° \n{-g 00 ) + f(t)yf^K*) (. . .) 

= j d i x^g-(f(t)+3Hf(t)+O(Sg 00 ,5K£))(...), (25) 

we can rewrite it into the second and higher order terms in cr, Sg 00 , and 8K p . Similar discussions hold for more general cases and we 
conclude that _F( 2 ) starts with quadratic terms of 8g 00 , cr, 5K pv , and 8R pupa . See also appendix B in |27|| . 
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where the dots stand for terms of higher order in the fluctuations or with more derivatives. When we rewrite the 
action in the unitary gauge in terms of the Goldstone boson 7r, the terms displayed in (|27[) are described by tt and its 
first order derivatives. In this paper, we consider the action up to the same order in derivatives of n and a. 

(2) 

Let us first construct the second order action. The second order action Sa containing a and its first order derivative 
can be written generally as 



S« = / d'x^g 



ai (*) iwa a , a 2(i) /a0 ^2 a 3(*) 2 , ,.\ qO 



(28) 



where we note that terms such as a(d ) 2 <j can be absorbed into other terms up to higher order fluctuations by 

integrating by parts. Here the second term leads to a non-trivial sound speed c a of a given by c£ = cei/{ai + a 2 )- 

(2) 

The second order mixing Si? is generally given by 



(2) 



s 



It is convenient to note the relation 



P^Sg 00 * + p 2 (t)5g 00 d°a + fo(t)6K£<T 



d±x^-gf{t)5K»{. ..)= I d±x^—g[f{t)d%. . .) (/(t) + 3Hf(tj) (...) + ^6g°°(. . .) + ^Sg 00 d°(. . .) + 



(29) 



/(*) 



(30) 



which can be obtained using (|25]l in footnote[6]twice. Here the last dots stand for higher order terms in the fluctuations, 
which can be written using 5g 00 , d°Sg 00 1 and SK^. Using this relation (|30|) . S^ x in (|29]l is rewritten as 

S£L = I d 4 x^\^(t)Sg 00 a + p 2 (t)5g 00 d°a + f3 3 (t)d°a - (&(<) + 3H(3 3 (t))a 



(31) 

where $\ = /3 1 +/3 3 /2, /3 2 = /3 2 +p 3 /2, and we dropped higher order terms constructed from a, Sg Q0 , d°5g 00 , and SKfi. 
In the following, we employ Eq. ()31l) as a definition of the second order mixing action S^ x . 

(3) 

The third order action S& of a is generally given by 

Si 3) = [ d 4 x^g-[ 7l (t)a 3 + 1 2(t)<J 2 d°<j + l3 (t)a(d°<j) 2 + l4 (t)(d°a) 3 + l5 (t) ad^a + l6 (t)d°ad^a\ , (32) 

(3) 

and the third order mixing S^( x is given by 

S ( l = I <?xV=gh {t)6g 00 a 2 + l 2 (t)Sg 00 ad° ( r + j 3 (t)S g 00 (d a) 2 +j i (t)d 6g™ad°a 



+ 75 (t) Sg 00 d,ad^ + %(t)(6g 00 ) 2 * + -f 7 (t)(Sg 00 ) 2 d°a 



(33) 



Here, it may be wondered if the terms such as 5K^a 2 , 6K^ad a, 6R 00 a 2 , and SK£6g ao a can appear at the same 
order in derivatives. However, they can be absorbed into other third order terms in (|3"3")l and the second order terms 
in (|28|) and (|31|) by integrating by parts as we did to rewrite (|29|) into (f33~|) . The term proportional to SRa 2 can also 
appear but such a term vanishes in the decoupling limit, so we do not consider it here for simplicity. 

To summarize, the most generic action in the unitary gauge can be written up to the third order fluctuations as 
follows: 



O _ a , c(2) , c(2) , o(3) , o(3) 
a — ag r av TO,, T J m i x < > J a > ^mix ! 

where 5 grav , S^ , S^l, S {3 \ and fl£^ are defined in (JUJ), (gSJ, USD, 422), and ([33]). 



(34) 



In the last of this subsection, we discuss the normalization of a. Using the field redefinition of a, we can drop some 
terms in the action and consider the following three normalizations. 
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(2) 

a. Normalization 1 : a\ + a 2 — 1- Let us first consider the kinetic term of a. The second order action S& of a 
can be expanded up to the second order fluctuations as 



d 4 xa 3 



( 21 

2 



o 2 -c 2 



(did) 1 a 3 — 3Hat4 — d 4 



Q'l + a 2 



(35) 



where the normalization factor a a is defined as a 2 = a± + a 2 . Although the factor a a is time-dependent in general, 
it can be taken unity by redefining a as a = a a a. Since the derivative of a can be written as 



n -It - cO - 

d A ,.<T = a a d^a - 6 —~o . 

at 



(36) 



the action still takes the form (|28| after the redefinition. 

b. Normalization 2 : = /J3 = 72 = 0. The action S in Eq. (|3"4"|) contains the following terms: 

f 



d^Xy 



-ff(t, aWd^adva + h(t, a)f 2 (t, a)d u a 



(37) 



where f±(t,a) and f2(t,&) are functions of t and a but do not contain derivatives of a. They are expanded in a as 
fi = a.\ — 275(7 + 0(cr 2 ) and fif 2 = /?3 + oyer + 72c 2 + 0(a 3 ). The function f 2 can be set to zero by redefining a 
such that 



f 1 {t,a)d ll a - S° f 2 (t,a) = h{t,a)d^a, 



(38) 



where h(t, a) is some function of t and <r but do not depend on derivatives of a . After the redefinition, the action still 
takes the form Q34p and the second order action containing a simplifies as 



Si? 



a l _a _ , "2(t)^0^2 "sM 2 



2 



p 1 (t)5g 00 a + ^ 2 {t)Sg m d a 



(39) 



c. Normalization 3 : a.\ = 1, 75 = 0. The action S in Eq. (|34[) contains the following term: 

f 



<i X\f—l 



--f 2 (t,a)g^d^d u a 



(40) 



where /(£, cr) is a function of i and a, and does not contain derivatives of a. The function is expanded in a as 
/ 2 (i,er) — ai(t) — 2j 5 (t)a + 0(<r 2 ). When a 2 is small enough, in other words, the sound speed of a is almost unity, 
it would be convenient to take f(t,a) to be unity. Let us define a as a = F(cr,t) such that dF/da = f. Since the 
derivative of a is given by 



9^ = 18^ + 6°^, 



(41) 



the action can be rewritten as 



d x\f— ( 



8F 



1 00 IdF 



at 



(42) 



which still takes the form (|34|) . Here F should be regarded as a function of u and t. It also should be noticed that, 
though we can set the function / to be unity, the terms proportional to d°a and g 00 appear. 



B. Action for the Goldstone boson and the decoupling limit 



In this subsection we construct the action for the Goldstone boson 7r and discuss its decoupling regime. As in the 
last section, we perform the time diffcomorphism (|10[) . Practically, it is realized by the following replacements: 



g oo ^ g oo + 2d o n + Q^d** , d°a -> d°a + d^n d^a , 
a ->cr, /(*)—>/(* + 7r) , J d 4 x\f~^g^ J d 4 x^J~^g. 



(43) 
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(2) (2) 

With these replacements, 5 gr av! S„ , and SV ; X are rewritten as 



Sgrav — J d X\J ~g 



M^R + M^H(t + tt) (g w + 2d% + d^d^ir) - M% x ( 3H 2 (t + tt) + £T(f + tt) 



j(2) 



J d 4 xy/^g 



M 2 4 (t + 7r) 
2! 

Q!l(t + 7r) 



(5g 00 + 20°7T + ^ttS'V) + 



3! 



(44) 



a 2 (t + n) , ou„\2 a 3 (t + 7r) o 



V + a4(t + tt)ct(5V + <9V) 

(45) 



/3i(t + tt) (c5<? 00 + 20°* + dp-icWir) a + p 2 (t + n) (Sg 00 + 2d°n + d^d^) (d°a + d^d^a) 
+ /3 3 (t + Tr)(d° (7 + 8^ a) - ($ 3 (t + tt) + 3H(t + tt)/3 3 (< + 7r))cr 



(46) 



(3) (3) 

The third order actions and S'j nix can be obtained in a similar way. 

In order to discuss the decoupling regime of the action, we first clarify in which regime graviton fluctuations become 
irrelevant to tree-level three point functions of tt. For this purpose, let us take the spatially flat gauge (fT5|) and use 
the ADM decomposition: 



ds- 



-(N 2 - NN l )dt 2 + 2N l dx i dt + a 2 (e 7 )y dx i dx j with 7li = = . 



(47) 



Here and in what follows we use the spatial metric hij = a 2 (e 7 )ij and its inverse h %3 = a 2 (e 7 )jj to raise or lower 
the indices of N l . The inverse metric g^ v are written in terms of N, N l , and hij as 



9 N 2 ' 3 9 N2 , 9 it- N 2 



(48) 



In this gauge, there are no second order mixing terms of 7r and jij because 7^- has two spatial indices and is transverse- 
traceless. Then, the tensor fluctuation 7^ does not contribute to tree-level three point functions of tt. Therefore, 
possible contributions of graviton fluctuations come only from the auxiliary fields SN = N — 1 and N l . As discussed 
in (30I ] . it is sufficient for the calculation of three-point functions to solve the constraints up to first order. Expanding 
the actions (|44|) - (|46| up to the second order in tt, SN, and N l , 



s. 



(2) 



d A xa 3 



- M^(3H 2 + c- 2 H)SN 2 - 2Mp 1 HSNd i N i + M^N l d t djN j - Mp^N 1 d 2 N l 



m 2 x h ( ^ 2 (d^y 



3M^H 2 tt 2 + AI^{2c- 2 Htt - 6HHtt)SN + 2M^HN l diTT 



2 TJAji', 



d 4 x a 3 
d 4 x a 3 



al a 2 -c 2 



(cV)< 



"3 2 

— a — oiaOO 
2 



2/3i (SN - tt) a - 2f3 2 (SN - tt) & - SN ( - (3 3 & + (3 3 a + 3H/3 3 a ) + (3 3 N l d % <r 

' p 3 TTG — 3Hf3 3 TTa 



~\~ $3 [ TT& 

the constraints are solved up to first order as follows: 

p3 



5N = -—tt-- 

H 2M 2 l H 



a, N l = a~ 2 diip 
H 



with ^ = a 2 d- 2 [c- 2 ^(HTT + HTT) + ^—a ^— a 



2M 2 l H 



-2 H 

c * H a + a 



2M 2 X H 



(49) 
(50) 

(51) 



(52) 



Here the sound speed c 2 of tt are defined as c 2 = HMpJ(HM^ — 2M|). The factors c 2 = ai/(a\ + a 2 ) and 
a 2 = ai + a 2 are the sound speed and the normalization factor of a, respectively. Using the canonical normalization 



n c ~ Mp^-H^c- 1 **, a c ~a a a, SN C ~ M P1 5N , JV* ~ M Pl N l 



(53) 
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and redefining the coupling constants /3i, /?2, and /?3 correspondingly as 



yi " o^Mp^-ff) 1 /^ 1 ' ^ ~ a(r M P i(-ij)V2 



ft, Pt ~ 7 ^, ln fa: 



a, 



■M Pl (-H)V*' 



we rewrite the constraints (1521) as 
N< 



e 1/2 § ( - + ^F7r c + /3?(7 C - /3 2 c (d c - ^a) - i/3 3 c d c + ljla c [ - 2^a + (<£ \)IH \n 



(54) 



(55) 



(56) 



_ 2 # - e 

where we have defined e = — c — ^ and rj = — in analogy with usual slow-roll parameters e and 77. It is manifest 

that (5-/V c and iV* are suppressed by the parameter e 1 / 2 and contributions from 5N C and JV* become irrelevant in 
the limit e — > 0. 7 In this limit, tree-level three point functions of ir are determined by the following action in the 
decoupling limit: 



S <r 



:(2) 



d 4 xa 3 



d 4 xa 3 



AM 4 

7T 



3M 2 ,H 2 tt 2 - d t ( tt 2 7t + MlH^^K - 3M 2 HHtt 3 



(57) 



(did) 2 Q!3 — 3Hct4 — Q.4 



\ , 2/1 2 \ ( ■ ■ 2 • (dindi<j) 
-a + a a (l - c CT )l TTcr - 



/ v or 



(.. d t nd t a\ (a 2 ) . 2 (a 2 c 2 ) (d.;cr) 2 d 3 2 

C^CtI 7T<T t. H — 7TfJ — - — 7T ^ (7 7T — a47T(TfT 



(58) 



d 4 xa 3 



(-2ft + ft)™ + (2ft - ft)** + ft^f^ - 3Hfo™ - A (V ^ 



C7+ (-2/3l +/3 3 )7r7T(T 

+ 3ft* 2 <7 - - ft^a + (2ft - ft)™* + kA^f) 3 (ffft + ^ft) tt 2 * 



S^= / d 4 xa 



(cV) 2 . (cV) 2 



(71 + -H72 + -72) o- 3 + (73 - 7 5 )<7o- 2 + (-74 + 7 6 )* 3 + 75 cr— ^ 7 6 (T 



? (3) 



d 4 xa 3 



2 / \ 2 ( ) — 2 — 2 

— 27 1 tt(7 + 27 2 ^tct(t + 2(— 7 3 + 75)71"* — 27 4 7rCT(T — 27 5 7r V 47 6 tt <t — 47 7 7T * 



(59) 
(60) 

(61) 



It should be noticed that non-trivial cubic interactions appear generically when the sound speed c a of a is small, 
is non-zero, or mixing couplings ft and ft exist as well as the sound speed c n of ir is small. 

C. Examples 

In the last of this section, we clarify the relation between our approach and models in the literatures. For this 
purpose, we first discuss the original model of quasi-single field inflation [2l|, and then, we investigate the effects of 
heavy particles during inflation. In the last of this subsection, a class of two field models will be considered. 



~ To be precise, we need to assume that e is small enough to vanish when multiplied by other parameters in the calculation. For example, 



we need to assume that e-^- <C 1. 
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1. Original model discussed by Chen and Wang 

As was reviewed in section [Til the original model [2l| of quasi-single field inflation is described by the following 
matter action: 



Smarter = / d A X^~g[ ~ \{R + *) ~ ^ - V S r{9) - V( X ) 



(62) 



The homogeneous backgrounds are given by 9 = 9o(t) and \ = Xo (constant), which leads to the action in the unitary 
gauge 56 = 6 - 9 Q = 0, 



Srr 



d x\J—g 



rf 4 



\{R + cxfg m 9l - \g^d^d v cj - V ar (9 ) - V( X o + <r) 
\RHlg m {VM + V( X o)) - \g» v d»o%o l(v"( X o) 9 2 )a 2 



R9£ 5g vv <7 



3! 



(63) 



where R = R + xo, & = X ~ Xo an d we used the background equations of motion. Using the equations of motion (1201) , 
the action (|63[) can be written in terms of the Hubble parameter H as 



J d^x^—g^Hg™ - M^(3H 2 + H ) - ±g^d,ad u a - \(v"{ X o) + 



2Mg 1 ff 
i? 2 



J? " " 3! J? 2 
which corresponds to the following parameters in our framework, 



6g 00 a 2 + 0(6a 4 ) 



(64) 



ai = l, a s = V"(xo) + — , Pi 



2M^H 
R 



7i = -^V"'(xo), 7i 



— (others) = 0. 65 
it- 2 



2. Effects of heavy particles 

Recently, it is argued that the existence of heavy particles can cause a non-trivial sound speed of effective single 
field inflation [13l420| . As was mentioned earlier, our framework is also applicable for such inflation models with heavy 
particles. In the following, we briefly show that it gives a simple explanation for the effective sound speed. 

As a simplest example, let us consider the following quadratic action: 



S 



d i x^/—g 



1 • .1 m 2 

-M 2 X R + M^Hg 00 - AI 2 ,(3H 2 + H) - -g^d.a d v a - —a 2 + 5g m a 



(66) 



Here m is the mass of a and /3 is the mixing coupling between the adiabatic mode and the massive particle a. We 
assume that time-dependence of m and (3 is negligible, and the mass of a is much larger than the Hubble scale during 
inflation, m 3> H. In such a regime, the kinetic term of a becomes irrelevant and the dynamics is determined by 



S 



d i xy/—g 



l -M 2 x R + M 2 x Hg™ M 2 } (3H 2 + H) - ^ 



P 



Sg° 



JUL. 
2m 2 



(Sg c 



(67) 



which implies that the perturbation a quickly responds to the variation of the adiabatic mode 5g 00 . Integrating out 
the massive particle <r, we obtain the following effective action for single-field inflation: 



S 



d 4 x 



hl 2 ,R + Al 2 ,Hg m - M 2 X {ZH 2 + H) + ^(Sg 00 ) 2 
2 2m z 



(68) 



As we discussed earlier, the last term gives the following non-trivial sound speed: 



-HM 2 X 



2P 2 /r 



(69) 
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3. A class of two-field models 



Let us then consider a class of two-field models described by the following matter action: 

d A x^\ - \lab{r)9^d^ a d^ b - V(</> a )] (a = 1,2), 



(70) 



where "f a b{4' a ) is the metric on the field space. This class of models were carefully studied in [31| and recently discussed 
in [TEj to investigate effects of massive particles during inflation. Suppose that the trajectory of the homogeneous 
background fields 0g (t ) is on a curve </> a = <fi a (X), and the background fields are given by 0g (i) = <fi(Xo(t)). We also 
assume that A > 0. Defining the coordinates (A, a) of the field space such that the curve a — coincides with the 
trajectory curve, the fields X(x) and <j{x) describe the adiabatic mode and the isocurvature mode, respectively. There 
are still many choices of the coordinates or degrees of freedom of the field redefinition. In the following, we consider 
two types of basis of the fields and discuss their properties in our framework. 

a. Orthogonal basis We first consider the orthogonal-basis. We can always take the coordinate (A, a) such that 

7Acr = 7o-A = . (71) 

By the field redefinition of A, we further require 7aa(A, <t = 0) = 1. In this basis, the matter action (|70[) is given by 



Smattcr — / d 



%V -9 



\lxx(X, a^c^AcU - i 7ffff (A, a)g^d^d v a - V(X, a) 



(72) 



and it can be expanded in the unitary gauge SX = A — Aq = up to the third order fluctuations as 



d x\f—g 



1\2 00 1 '\2( \l r 00 i 1 \2/ \//_ 2 M2/ \ii r 00 2 . 1 {2/ 3 

- 9 V9 - 2 A o(7aa) <5s cr + -A o ( 7 aa) o 0- - ^X {-y\\) dg a + — A ( 7 aa) cr 
- ^(7 w )o^V^ ff " \<a<,„)'<i<r9* v d^du<r -V - \v "a 2 - ±V "'a 3 



(73) 



where we have used the equation of motion for a. In this subsection, we write derivatives of the metric and the potential 
evaluated at the classical value, for example, as (7aa)o — <9 ct 7aa(A, ct)|a=Ao,<x=o an d V Q ' — d a V(X,o-)\\=\ 0t(T =o- Using 
the equations of motion for graviton, 



1 ■ 



A^ , M^(3H' +H)=V , 



(74) 



we can rewrite (1731) as 



d xyj—g 



M 2 l9 w H - MUSH* + H)- ^( 7 ™)o S^crcU - ^( U " + M^H^xx)' ' W 



1 



1 



+ M^H( 7xx )' 6g 00 a - -( 7 ™) a g^d^d v a - - V "' + M^H{lxx)o °* + -M^if ( 7 aa)oVV , (75) 



1 



which is described in our framework as 



a\ = (7<t<t)o , a 3 = Vq + M| 1 F(7aa)o , Pi = m piH(ix\)'o ■ 



7i 



G 



AAJq 



75 = -^(7««r)o . 7i = 2-Mp 1 ^(7 AA )o , ( others ) = . 



We notice that this basis corresponds to the second normalization of a in section IIII Al 
b. Canonical a basis We then choose the coordinate (A, a) such that 

laa{X,a) = 1 , 7aa(A, a = 0) = 1 , 



(76) 



(77) 



where the normalization of a is always canonical and that of A is canonical only on the trajectory curve. In this basis, 
the matter action ([70)1 is given by 



= I d x\J~^g 



\lXx{X, G)g^d^Xd v X - 7A(7 (A, <j)g^d^Xd u <7 - \g^d^d v a - V(X, a) 



(78) 
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and it can be expanded in the unitary gauge 5X = X — Xq = up to the third order fluctuations as 



5 n 



d xy/—g 



1 \2 „00 , 1 \2 



"~3 



Ao(7A.)odV - Ao( 7 A CT )oaa°a - iAo( 7 A CT )oV 2 a°a - ^d„ad v a 

T/ T/7_ ^T/"„2 ^T/"'^-3 



(79) 



where the equation of motion for a implies Vq = ^Xq(j\\)' — ( 3H Ao(7a<t)o + Ao(7 A(T )o + Ao 9t(7Acr)o ) • In terms of 
the Hubble parameter H, we can rewrite it as 



= / d 4 xy/^g 



M 2 ,g m H - M^(3H 2 + H) - ^M^-H) 1 ' 2 ^)^* - (V% + M^H(~/ xx )' )a 

- \g^d^d„o - i(y " + M| 1 ij( 7 AA) , o')^ 2 - x/2M P1 (-H)^ 2 ( 1Xa )' ol jd a a + M^H{ lxx )' Q 8g m cj 



which is described in our framework as 

ai = l, a 3 = K' + m pi^(7aa)o, a 4 = -^A/ P1 (-ij) 1 / 2 ( 7Aff )J ) , ft = M^H(-f XX )' , 

& = -a/2Mpi(-F) 1 / 2 ( 7 a ct )o , 7i = ~\ (C + M P1 i/( 7 AA)£') , 72 = -^pi(-# ) 1/2 ( 7 a^' 
1 



(80) 



7i = ^£(7^ , ( others ) = . 



(81) 



The above action apparently includes terms linear in er. However, it can be easily shown that such terms turn out to 
be more than second order after integration by parts. We notice that this basis corresponds to the third normalization 
of a in section IIII Al 



IV. POWER SPECTRUM 



In this section we calculate the power spectrum for a class of quasi-single field inflation models in our framework. 
In the following, we take the decoupling limit e — > and use the action for the Goldstone boson tt. It is also assumed 
that the background trajectory satisfies eCl. Up to the second order fluctuations, the action for 7r and a takes the 
following form in the decoupling limit: 



S 



d 4 xa 3 



M P1 H 1^2 _ c 2 



2 (^ CT ) 2 2\ , 5 ■ , o ■ ■ , 5 2 di'*d l CT 



mX + /5i7TCT + ft Tier + ft<- 



(82) 



where we have dropped sub- leading terms ~ Mp^Hir 2 and ~ Hirer in the regime e -C 1, and have defined 
m , = a s -3Ha 4 -a^ p 1 = _ 2 !3 1 +$ 3 , ft = 2ft ft , ft = c~ 2 fo . 



(83) 



In this paper we concentrate on the case where the mixing couplings can be treated as perturbations although they 
can give large contributions in general. As pointed out in [211 ] . the perturbativity of mixing couplings is justified even 
in the case when they are large only in a sufficiently short period of time, which is realized, for example, in the case 
of the sharp turning background trajectory [l4l - fl7l [l9l[20j. 

Let us then calculate the power spectrum using the in-in formalism. Before going to concrete models, we introduce 
a general expression for the power spectrum up to the leading order corrections of the mixing couplings. As in the 
usual case, the fields ir and a are expanded in the Fourier space as 



TTk = u k a k + u* k at k , cr k = v k b k + vt bK 



(84) 
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with the standard commutation relations 

[<z k , 4,] = (27t) 3 <^ (k - k') , [6k, bl] = (2irf6& (k - k') . 
Here the mode functions u k and v k satisfy the equations of motion in the free theory and depend on k = |k| 



ii k + H(3 + 2e + fj)u k + ^—u k = 0, v k + H(3 + 2 S a )v k + m* + )v k = with 5 a = 

Their normalization follows from 

2Mp l H 2 ea 3 (u k u* k - u k u%) = i, a 2 a 3 (v k v* k - v k v%) = i . 
Using these expressions in the interaction picture, the expectation value of 7r k (i)7r k ' (t) is given by 



Ha„ 



(85) 



(86) 



(87) 



M*)*v(*)> 



Texp(ijf dt'H^t')} 



7Tk(i)7r k /(i) 



Texpf -i / dt'HAt' 



10) 



to 



(o|7r k (i)^(t)|o)+ / / <Mo|#&(*iM*)*v W^i7 x (<i)|o) 



2Re 



dh / ^2 (0| 7T k )7T k , ) i?i 2 iL (* 1 ) |0> 



(89) 



(2) 

where the dots stand for the higher order terms in the couplings and * s ^ ne secoil d order interaction terms 



d 3 k r~ ~ k 2 1 

a 3 (t) /3i7T k cr_ k + /3 2 7r k (T_ k + /3 3 c^^-7r k CT_ k (i). 



(27T) 



(90) 



In terms of the mode functions and the couplings, the general form of the two point function (|89j) is given up to the 
leading order corrections of the mixing couplings by 



(7T k 



(*)**(«)) = (27T)V 3 )(k + k'K(t)u fc (t) 



1 + Ci + c 2 



where Ci and C 2 are defined by 



d=2 



dti a 3 



to 



Pi U k V k + p 2 U k V k + /?3C^-2 U kVk 



(tl) 



(91) 



(92) 



Co = -4Rc 



u|(t) 
l«*(*)l s 



/3i u^wfe + p 2 u* k v k + p 3 c^^u* k v k 



(ti) 



dt 2 a 3 



(*2) 



(93) 



It is convenient to rewrite C = Ci + C 2 as follows: 



C = 4Re 



dh a 3 



(~ ~ k 2 

( /3i Ufcffc + p 2 ii k v k + P-icl^r u k v k 



uj{t) 

\Mt)\'< 



(~ ~ k 2 

( p 1 ii* k v k + p 2 u* k v k + p 3 cl-^u* k v k 



(tl 



dt 2 a 3 



k 2 1 

px u* k v* k + p 2 ii k v k + focl-^ u* k v* k \ (tj,) 



(94) 



Since the scalar perturbation C, is given at the linear order by £ = —Hit, wc obtain the expectation value of Ck(i)Ck' (i) as 



(Ck(*)Ck'(f)) = (27r) 3 ^ 3 )(k + k')^P C (fc) : 



where the power spectrum 'Pf(fc) is given by 

V c (k) 



H 2 (t)k 3 
2ir 2 



■«l(t)u k (t)(l+C) 



(95) 



(96) 
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The factor C can be considered as a deviation factor from single field inflation. Here it should be noticed that in 
the derivation of the above general expression we assumed only e <C 1, e -C 1, and the perturbativity of the mixing 
couplings. In principle, we can calculate the power spectrum for any models satisfying those three conditions. In 
the rest of this section, we first perform the calculation concretely for the case when the time-dependence of mixing 
couplings is irrelevant (constant turning trajectory), and then we discuss the qualitative features of the case when the 
mixing couplings are large in a sufficiently short period of time (sharp turning trajectory). 



A. Constant turning trajectory 



In this subsection, the power spectrum is calculated in the case that the time-dependence of mixing couplings is 
irrelevant. To make the calculation tractable, we assume that the time-dependence of e, e, a a , c a , and m a is negligible 
and we use the de-Sitter approximation. In this approximation, the equations of motion (|86p for the mode functions 
Uk and Vk can be written as 



9 9 m 2 

II I i 2 i 2 n II I i 2 ; 2 i (j n 

u k - ~u k +c n k u k = 0, v k - -v k + c a k v k + -rr^vk = , 



(97) 



where the conformal time dr = a 1 dt is given by t = —l/(aH) in the de-Sitter approximation and the primes denote 
derivatives with respect to r. The equations (!§?]) can be solved as follows: 



Uk 



2A/pie 1 /2( CT fc)3/: 



■(1 + ic^kr^e 



—iCtt kr 



1 



2A/pi? 1 /2( C7r / c )3/: 



■(l-ix)e ix , 



2a a 2a (J {Cirk) 4 i 2 



(98) 
(99) 



where /c n , and we chose the Bunch-Davies vacuum for 7r and a. The function H„ = J v + iY v is 

the Hankel function and v is defined as 



rn a <-H, v = l ^ — - - for m a > -H . 



The time derivatives of u k and Vk are given by 



Uk = —Htu'u 



H 



Vk = —HtvI = ie 



2M P1 ? 1 /2(c 7r fc)3/2 



,x 3 / 2 ((3/2 - u)H^{r s x) + {r s x)H^\{r s x)) , 



where we used the identity zd z H„ = zH^\ — vH^\ Therefore, the factor C defined in takes the form 



(100) 

(101) 
(102) 



C = e^ {v - u ' ) 



x Re 



4o£ M^i-H) 



dx\ 



+ {v - 3/2)*/3 2 - ft) (M{x)* - Ai(xi)) - /3 2 (a 2 (x)* - A 2 ( Xl )) + &(a 3 (x)* - l 3 (zi)) 



x / dxi 



+ {y - 3/2)^ 2 - ~^A x {x 2 ) - p 2 A 2 (x 2 ) + PMx 2 ) 



(103) 



where we set to = — oo and t = oo, and we defined 

A 1 (x) = x-^^H^inx) , A 2 (x) = r.x^e^H^r.x) , A 3 (x) = ^"^(m) , 
A! (a:) - x~ 1 ' 2 e ix H£\r B x) , A 2 (z) = r^/V^l^z) , A 3 (x) = ix^e^H^ (r s x) . 



(104) 
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When the time-dependence of mixing couplings is irrelevant, mixing couplings /3's can be evaluated at the time of 
horizon-crossing r = — (c^fe) -1 . In such a case, C is given by 



Re 



ft/J2 + {y- 3/2)/3 2 - ft In + 0i/H + (y* - 3/2)ft - ft)(-ft 2 12 + ft2i 3 ) 
- ft ((ft/# + 3/2)ft - ft)2 2X - ft 2 22 + ft 2 23 ) 
+ ft ((ft/J2 + (y - 3/2)ft - ft)2 31 - ft2 32 + ft2 33 , 




ft 2 ft 2 



Pi h 



Pi C 33 + g ft Cia + g ft C13 + ft ft C 23 



where 1^ is integral defined as 



-e 2 



7r(f — 1/*) 



dxi 



\Ai(xi) 



Ai{ Xl ) 



dx 2 Aj(x2) ■ 



(105) 
(106) 

(107) 



and Cij is given by 



C n = Re 



X11 



C 2 



Re 



W-3/2\ 2 l n 



3/2)2i 2 -(z,-3/2)2 21 +2 22 



C 33 = Re[Jii+2'33-Ii3-X3i] , Ci 2 = Rc[(^ + ^-3)Iii-Ii2-I 2 i] , C13 = Re[-22n+2i 3 +2 3 i] , 
C23 = Rc [-{v + v* - 3) 2n + 2 12 +2 21 + {v* - 3/2) 2 X3 + (1/ - 3/2) 2 31 -2 23 -2 32 ] . (108) 

Here the explicit form of the power spectrum is given by 

4 (M 



V c (k) 



H 2 
8Tr 2 M^€c n 



1 



Cn + ft 2 C22 + ft 2 C33 + % ft C12 



if 



1 ft C13 + ft ft C 2 3 



(109) 



where H, e, c^, q ct , ft, and r s are evaluated at the time of horizon-crossing. Then, the calculation of the power 
spectrum reduces to the evaluation of Re[2^ + 2 3 -,], Im[2i 2 — 2 2 i], and Im[2i3 —231]. 

As is understood from the definition, 2jj's and C^-'s are functions of m a and r s = c a / 'c„. For general value of r s , it 
is difficult to perform the integrals analytically and we performed numerical calculations by contour deformations (see 
appendix |A1 for details). For the special case r s = 1, however, it is possible to perform the integrals 2^'s analytically 
by extending the results in [22I ] , and the results are summarized in appendix [Bj In such a way, Cij 's are calculated 
and the obtained results are summarized in figure [T] and figure [2J We find that they monotonically decrease in m a 
for fixed r s , but they are not monotonic for fixed m a . It is also notable that C 22 does not vanish even in the limit 
m a — > 00 contrary to the usual expectation that heavy particles decouple for finite mixing couplings. In the last of 
this subsection, we give an explanation of the peculiar features of heavy particles. 

In the heavy mass limit 3> H and m a ^> a _1 c cr fc, the mode function Vk can be expressed as 



1 



Vk 



a -3/2 e -im a t i 



(110) 



where the normalization of vu follows from (|87[) . Then, the mixing couplings ft and ft are suppressed by the 
prefactor 1/ \jm a and these mixings do not contribute to the deviation factor C in the heavy mass limit. On the other 
hand, in the case of the ft mixing coupling, such a suppression does not occur because ik takes the form 8 



Vk 



V2a a 



a -3/2 e - lm „t 



(111) 



Apparently, the expression seems to imply a divergence in the mixing. However, it turns out that such a divergence does not 

occur because of the integration over time. 
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Cn 




m a /H 



m a /H 



C 



:-!:! 




m a /H 



1 i i 8 f f » 




m a /H 




rn a /H 



m a /H 



FIG. 1: Cij's for fixed r s = Ccr/c^. The dots are numerical results for r s = 0.1 (red), 0.3 (orange), 1 (yellow), 3 (green), and 10 
(blue) . The curves are analytic results for r s = 1 obtained in the next subsection. 

More explicitly, it is possible to evaluate C22 in the heavy mass limit. In this limit, Vk and Vk are related to each other 
as 



and C22 is related to Cn as 



v k = -im a Vk , 



r - m °r 
~ IP 1 



As was discussed in [THHHi Cn for r s = 1 in the heavy mass limit is analytically evaluated as 9 

1 H 2 



Cn 



4m2 



(112) 



(113) 



(114) 



We then obtain 



C22 = t in the heavy mass limit. 



(115) 



9 The factor C(fi) in [Hi |22f corresponds to our Cn for r s = 1 and [i ~ — — in the heavy mass limit. 
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FIG. 2: Cy's for fixed m a . The dots are numerical results for m&/H = 1.2 (red), 1.6 (yellow), 2.0 (green), and 4.0 (blue). 



Note that, in the heavy mass limit, Vk and therefore C^-'s do not depend on r s . The result (|115[) is consistent with 
our numerical and analytic results and we conclude that the 02 mixing coupling contributes to the deviation factor C 
even in the heavy mass limit. It would be important to understand this feature more physically. 



B. Qualitative features of sharp turning trajectory 

In this subsection, we discuss qualitative features of the case when the mixing couplings are large in a sufficiently 
short period of time (sharp turning trajectory). As a simplest example, let us first consider the case when the /?i 
coupling is the only relevant mixing coupling and it is proportional to a delta function: 

h =P.5(t-t.) =/3*a- 1 (Q5(T-T*), /3 2 = /3 3 = 0, (116) 

where t* is the time of sharp turning, r* is the corresponding conformal time, and the mode /c* crossing the horizon 
at t = t* is given by k* = —1/(c v t*). We assume that /3* is small enough to be treated perturbatively. In this case, 
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the deviation of the power spectrum from that of single field inflation can be written as 10 

C = 2fc Re [a 6 (u* k - u k ) u k v k v* k (U)] , (117) 

which is related to the power spectrum by (|96p . For simplicity, suppose that time-dependence of m ai a a , c a , and c w 
is negligible at the time of sharp turning and the mode functions u k (t t ) and v k (t*) are given by 

(118) 

v k (Q = -ie*™+i* ^ H xl^HpXrsX.) , (119) 
where = fc/fc* and parameters such as m a are evaluated at the time of sharp turning. Then, (|117[) is given by 

2 

C = 0l * 9 C " -M x,) with F v {x^ = ei^-^x,sin 2 xAH^\r s x,)\ 2 . (120) 
For k <C fc* or s, < 1, the asymptotic behavior of J- V {x*) is 

^T{vfr- 2v xl- 2v for m<T (U) < |iJ , 

■7v(*.)~< ^ 2 (121) 

- e^ 7rl/ 2 l T(i/)(r s a; >t )~" + e-* 7r,y 2- l/ r(-i/)(r s a;,) l/ a* for m CT (t»)>-#, 



2 



which is consistent to the intuition that modes outside the horizon at t = i* are not much affected by the sharp 
turning. For k 3> fc» or x* 3> 1, it reduces to 

2 1 

J"„(a:*) ~ — sin 2 = — r ; T 1 (l — cos 2a;*) . (122) 

7T 7T 

This kind of oscillating behavior was also found in (TB-E3, Ell ■ The turning trajectory generically oscillates around the 
turning point and the mixing couplings at the turning point become more regular than delta functions. In such a case, 
it is expected that the oscillating behavior of short modes c^k > — 1/r* begins damping at some scale characterized 
by the oscillation of the trajectory. Let us next confirm such a behavior explicitly for a concrete example with a finite 
width of turning. We consider the following 0i profile: 

Pi = < 2 2 (123) 

otherwise , 

where we normalized /?* so that it reproduces the coupling (|116p in the limit At — > 0. For this class of couplings, we 
define T v by 

'-X^lds)^- (124) 

Since it is difficult to calculate J- v analytically, we can confirm the expected damping behavior of T v by numerical 
calculations and the results are given in figure [3] To summarize, C vanishes in the long mode limit c^k <C — 1/t„, 
oscillates for short modes c^k > — l/r», and damps at some scale characterized by the oscillating trajectory around 
the turning point. 

It is straightforward to extend the above discussions to the case with non- vanishing /3 2 and p 3 . We generically 
expect to find a similar behavior of C: vanishing in the long mode limit c^k -C — 1/r*, oscillating for short modes 
c^k > —1/r*, and damping at some scale characterized by the oscillating trajectory around the turning point. 



In our calculation wc have treated the mixing coupling /3» as an interaction. In such an interaction picture, it is manifest that the 
deviation from single field inflation starts from the second order in /3». On the other hand, it is also possible to treat the mixing coupling 
as a part of the kinetic and mass terms. In that picture, the commutation relations 185B of creation and annihilation operators are 
affected by the mixing as well as the mode functions u k and v k are modified. However, in some literatures, these modifications are not 
taken into account adequately and the deviation from single field inflation is calculated to start from the first order in /3* . 
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FIG. 3: Ti> for fixed c x i:*Ar and = c CT 
. The left figure is for m CT = H and the right one is for m a = 3-ff. The dots are numerical results for c„fcAr = 5 (red), 10 
(green), and 15 (blue). The curves are analytic results in the limit of c^k, A — > (delta function limit). 



V. THREE POINT FUNCTIONS IN THE SQUEEZED LIMIT 

In this section, we discuss the momentum dependence of three point functions in the squeezed limit. We take the 
decoupling limit and assume that time dependence of a's, /3's, 7's, and H is negligible. Under these assumption, the 
second order action is given in ([82"]l and ([55]). and three point vertices are given by 



5 (3) 



d xar 



(M| ljf j(c- 2 - 1) + - 4M 3 V 



)tt 3 + M^H(c~ 2 - 1)tt^- + {-Pi + 4%)n 2 <j + (3/3 2 - 477)71^ 



. diirdicr „ (diir) 2 „ (Ba-k) 2 . 9 „ , . . 

2/3 2 tt 9 + /3i^V-a - ^i-t-^-ff - 27 1 7tct 2 + (-a 4 + 2j 2 )w<t& 



or 

.2/1 „2> 



^(1 



273 +275 ™ -275* 



. (5icr) 2 9j7r9i<T 2 2 diTrdia . . 

; h «4 5 — cr - a CT (l - c CT ) 5 — a - 27 4 7rerer 



3 2 • / n • 2 (dicr) 2 , \ . q . (<9 ?; er) 2 
+ 7x0- - 7 2 (T a + (73 - 75)o"cr +750- 5 l-(-74 + 76jc - 76<7 o - 



(125) 



In the following, we discuss what kind of momentum dependence in the squeezed limit appears from the above three 
point vertices. 

As an example, let us start from the case when the mixing coupling $i"frcr and the three point vertex A'jQ^a are 
relevant. In this case, the three point function of 7t takes the form 



(7rki(*)7rk 2 (*)TTk 3 (*)> 

^ 76 (27T) 3 ^ 3 )(k 1+ k2+k3) 

3 M 3 1 e3/ 2 (c^ 1 ) 3 / 2 (c,fc 2 ) 3 / 2 (c^ 3 ) 3 / 2 



Re 



dtia u* kl {ti)u* k2 (ti) 



x |^fc 3 (*i) J dt 2 a 6 u k3 {t 2 )vk 3 (t2)~v* k3 (t 1 ) I dt 2 a d ul 3 (t2)v k3 (t2)-v k3 (t 1 ) / dt 2 a 3 u* k3 (t 2 )v* k3 (t 2 ) 
+ (2 permutations) . 



(126) 



Here there are three terms in the curly brackets, and they correspond to the Feynman diagrams in figure^! respectively 
In the squeezed limit, k\ = k 2 = k and k^/k = /s<l, the long mode k% crosses the horizon much earlier than the 
short modes k\ and k 2 . Since the interactions are considered to be relevant around the horizon, it is expected that 
the relevant contribution arises from (ti,T2) ~ (— 1/fci, — 1/^3), and therefore, the middle term in the curly brackets 
becomes irrelevant in the squeezed limit. In fact, we can confirm this expectation explicitly from the expression (|126|) . 
and the integrals in the curly brackets can be written at the leading order in k as 



dti a 3 ii* kl (t\)u* k2 (ti) 2i Im 



v* k3 (h 



dt 2 a 3 Uk 3 (t 2 )v k3 (t 2 ) 



(127) 
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FIG. 4: Feynman diagrams for the first term of (|126[1 . The solid lines denote the propagations of n and the dotted lines denote 
those of a. 



We notice that the t2-integral is ^-independent and A^-dependence of the total integral appears only via f/ji (ii) 
originated from the three-point vertex. Then, the momentum dependence of the first term in (|I26|) is given by 



the first term in (|I 261) 



K~ 3 / 2 ^ k- 



k z / 2 k 6 sinfzi/ log k + 8 V ] 



for m a < -H , 
lor m a > -H , 



(128) 



where 5 V is a i/-dependent phase factor. We note that the only information necessary to derive f|l 28[) was the fact that 
the field of momentum k% in the three point vertex takes the form a. For the other two permutation terms in (|126[) . 
the A^-dependence of the integral is determined by u* k3 originated from the three point vertex, and their momentum 
dependence is given by 



the other two permutation terms in (|126|) oc k k 



(129) 



Then, the first term dominates for small k. Therefore, when the mixing coupling j3\ and three point coupling 7 6 are 
relevant, the momentum dependence of scalar three point functions in the squeezed limit is given by 



, K-^ 2 -»k- 6 



for m a < -H , 

tor m a > —ti . 
2 



(130) 



It is straightforward to extend the above discussion for more general cases. First, for general Pi's, we can show 
that when the mixing couplings convert 7t of momentum k% = nk to a, the ^-integral becomes ^-independent in the 
limit K <C 1 and the three point vertex determines the K-dependence of the diagram. Second, the only information 
necessary to obtain the momentum dependence is the form of the field of momentum fc 3 in the vertex. In the previous 
examples, when it takes the form a in the vertex the diagram was proportional to KT i l 2 ^ v or k -3 / 2 sin[wlogft + S v ], 
and it was proportional to k" 1 for tt. More generally, we can obtain the relations in Table I between momentum 
dependence of the diagram and the form of the field of momentum k% in the three point vertex. Here it should be 
noted that the ^-dependent phase factor 5 U depends on the details of mixing couplings. Finally, as discussed in the 
previous example, momentum dependence of the contribution from each vertex is now identified for k -C 1, it is 
straightforward to obtain momentum dependence of the contribution to scalar three point functions from each three 
point vertex displayed in (|125[) . The results are summarized in Table II. Here note that although the contribution 

from the 7r — 1 —= — vertex seems to be proportional to K~ 2 k~ 6 apparently, explicit calculations show that this kind of 
leading contribution vanishes and the three point functions begin with terms proportional to K _1 fc -6 . 

As we have seen, the momentum dependence of scalar three point functions in the squeezed limit has robust 
information about mass of a and three point vertices. 



VI. SUMMARY AND DISCUSSION 



In this paper we discussed quasi-single field inflation based on the effective field theory approach. We first con- 
structed the most generic action in the unitary gauge based on the unbroken time-dependent spatial diffcomorphism, 
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form of the field of momentum k$ 


momentum dependence of the diagram 


7T 




7T, 7T 




din 
a 


K- 2 k- 6 


a, a 


1 


' K~ 3/2 - v k- 6 for m CT < 
K~ 3/,2 fc~ 6 sin[izv log k + S v \ for m a > 


dicr 
a 


! 


' K- 1/2 - v k- 6 for m CT < 
K _1/,2 fc -6 sin[i// log k + for m a > —J? 



TABLE I: Momentum dependence of the diagram 



three point vertices 


momentum dependence 


71 > n a 2 


AC-ifc" 6 


TV a, TV a, Tva , naa, ty<j , iraa, 

3 2. .2 (fta) 2 .3 .(<V) 2 
CT , CT (J, GO , CT , a , (J 

a z a z 




' K- 3/2 -"fc- 6 for m a < 
K~ 3/2 fc~ 6 sin[ii/logK + tv] for m a > 


. diTvdiCT 
W 2 




(din) 2 (dtn) 2 . dindia dindia . 

2 °"> 2 °"' 2 °"> 2 °" 

a- 3 a' 1 


1 


' K~ 3/2 -"k- 6 for m a < V2H 
K K- 2 k~ 6 for m a > V2H 


. (d t a) 2 

71 2 

a 2 


1 


' K- 1/2 -"k- 6 for m„ < V2H 
, K _1 fc- 6 for m CT > V2H 



TABLE II: Three point vertices and momentum dependence 



and then constructed the action for the Goldstone boson by the Stiickelberg method. Its decoupling regime was also 
discussed carefully, and the action in the decoupling regime implies that non-trivial cubic interactions generically ap- 
pear and large non-Gaussianities can arise when the sound speed c a of a is small, is non-zero, or mixing couplings 
/3i and /?2 exist as well as the sound speed c w of it is small. Using the obtained action, two classes of concrete models 
were discussed: the constant turning trajectory and the sharp turning trajectory. 

In the constant turning case, we first calculated the power spectrum of scalar perturbations numerically for general 
values of r s = c^jc^ and analytically for the special case r s = 1. It was found that when the tt<7 coupling is relevant, 
the massive scalar field a affects the dynamics of the Goldstone n even in the heavy mass limit m CT 3> H contrary 
to the usual expectation that heavy particles decouple for finite mixing couplings. We also discussed the momentum 
dependence of scalar three point functions in the squeezed limit for general settings of quasi-single field inflation. It 
was shown that the momentum dependence is determined only from the cubic interactions and the cubic interactions 
were classified into five classes. The three point functions in the squeezed limit take the intermediate shapes between 
local and equilateral types when the mixing couplings are relevant, and this kind of momentum dependence charac- 
terizes quasi-single field inflation. Recently in [32J , the detectability of such a momentum dependence was discussed 
for some cases. It would be interesting to discuss the detectability of the momentum dependence in the form of 
K -3/ 2 s in[ii/ log k + 5 U ], which arises in the second class with m a > %H. It is also important to calculate the full 
bi-spectrum for general settings of quasi-single field inflation. 
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In the sharp turning case, we made a qualitative discussion of the power spectrum. It was found that the deviation 
factor C from single field inflation vanishes for long modes c^k <C — 1/r*, oscillates for short modes c n k > — 1/r*, and 
damps at some scale characterized by the oscillating trajectory around the turning point. Since our framework makes 
the contributions from the mixing couplings clear, it would be useful to discuss more on the sharp turning trajectory. 

Our framework can be considered as a starting point for systematic discussions on multiple field models, and there 
would be a lot of applications such as those mentioned above. We hope to report our progress in these directions 
elsewhere. 
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Appendix A: Numerical calculations of power spectrum 

To perform the integral (|107p numerically, there arc two technical obstacles [U [22[ : spurious divergences at x = 
and oscillating behaviors at x = oo. 

The integral (|107p contains two integrals, J °° dx% Ai{x{)* dx2 Aj{x 2 ) and J °° dx\ Ai{x{) dx 2 Aj{x 2 ), and 
each integral diverges for some parameter region. However, we can show that such divergences cancel out in the 
calculation of C^ 's. For example, let us consider C\\ for real < v < 3/2. The asymptotic behavior of each integral 
in In around Xi = x 2 = is given by 

~ f dxi x- 1/2 - v [ dx 2 x- 1/2 - y ~ Q l - 2 \ (Al) 

Jo J x x 

which diverges for v > 1/2. We first notice that this kind of leading singularities cancel out between two terms and 
the asymptotic behavior of the total integral In is given by 

~i f dx lX \ ,2 - y f dx 2 x~ 1/2 - v ~i0 2 - 2 " (A2) 

JO J xi 

up to some real constant number. Although it is still singular for v > 1, such singular contribution is pure- imaginary 
and does not contribute to Re[2n]. The higher order terms are finite for v < 3/2 and hence we conclude that 
Cn = Re[2n] is finite. In a similar way, we can show that Rc[2jj + Zji], Im[Ii2 — ^21], and Im[2i3 — 131] are finite for 
< v < 3/2 and v = pure- imaginary, and therefore all the C^-'s are finite. To avoid this kind of spurious singularities 
in the numerical calculation, we introduce a IR cut off cir,: 

/>oo />oo 

Iii = -e^ v - v 'n dxdAiixi)* - Mxi)) dx 2 A j (x 2 ), (A3) 

Jem K 7 J x x 

where we set £ir = 10~ 10 in our calculation. 

As is usual in the Feynman diagram calculation in the momentum space, the integral lij oscillates at x = 00 
because of the oscillating behavior of the mode functions Uk and Vk, which makes the numerical calculation difficult. 
Following [ll|, [22|, we perform contour deformations to avoid this kind of technical difficulties. Let us first consider 
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the integral Rc[lij It is convenient to rewrite it as follows: 

" />oc />oo />oo />oo 

Re[l ij +l ji ] = -e^ u - u *)Re / d Xl A^)* dx 2 A j (x 2 )+ dxiAjfa) dx 2 A l {x 2 f 

.Jem J x\ Jem "Xi 

i pOC pOO poo 

dx\Ai(xi) j dx 2 Aj(x 2 ) - / dxiAj(xi) / dx 2 A t (x 2 ) 



Re 



tin 

oo 



oo 



dxA l (x)\ / <fyA,-(y) 

/•oo 

dxi^(xi) / dx 2 Aj(x 2 ) - / dxiA,-(xi) / dx 2 Ai(x 2 ) 



cm •'ai Jem J xx 

The first term in the bracket can be Wick rotated without crossing any poles as 

J dx Ai (eiR + ix) J J dy Aj (e m + iy) , 
and the last two terms are Wick rotated as 

/>OC pOO />OG 

dxAi(e m +ix) / dy Aj (e m + ix + iy) + / dx Aj(e m + ix) / dy Ai(ei K + ix + iy) , 
Jo Jo Jo 

Then, we obtain the following expression of Rc[Xij 



(A4) 



(A5) 



(A6) 



Re[2y +%] = -e^-^Re 



dx(Ai(em + ix)*+Ai(em + ix)) / 4?'( e iR + to + iy) 



+ (i o j) . (A7) 



To avoid the singular behavior around a; = discussed above, we further modify the contour as follows: 



i dx[ Ai(e m + x)* + Ai(e m + x) ) / dy A 3 (e m + ix + iy) 



+ J dx(Ai(e m + l + ix)* +Ai(e m + l + ix)j J dy Aj(e m + 1 + ix + iy) 
By performing similar contour deformations, Im[Iy — Xji] can be also expressed as follows: 



+ i) ■ (A8) 



Im[Iy - = -e 



+ -e^-^)lm 
4 



dx(A:(em + a;)* - ^(eiR + .T)) / dy Aj (e m + ix + iy) 
K ' Jo 

dx( Ai(ejn + 1 + ix)* + Ai(em + 1 + ix) J / dy A,- (e m + 1 + ia; + iy) 



- (i <-> j) • 

The expressions (|A8[) and (| A9|) are used in our numerical calculations of the power spectrum. 



(A9) 



Appendix B: Analytical calculation of power spectrum for = c a 

In this appendix we calculate the power spectrum for the case r, = c^/c^ = 1. For this class of models, we can 
analytically calculate the integrals I%j 's by extending the results in |22| ■ We first introduce a function A(£, v, x) defined 
by ' 



A(l,v,x) =x- 1 * +t e lx Hl 1 \x) 



(Bl) 



In terms of A, Ai can be written as 



At(x) = A(0,v,x) , A 2 (x) = A{l,v - l,x) , A 3 (x) = iA{l, v,x) 



(B2) 



24 



and hence the x-integrals in the first term and the ^-integral in the second term of (|107|) reduce to that of A. 
Similarly to the case in [22| , the indefinite integral T> of A can be expressed using hyper-geometric functions as 



V(£,vx) 



2 u x^ +l ~ v T(v) /l 1 3 \ 

■ ,t^ f ' 2F2 [--v,-+l-v;-+l-v,l-2v;+2ix) 
111(2 + £~ v) \2 I I / 

, ^^x^+Ti-v) (I _ 1 , t , 3 
-2^2 



(- + v,-+£ + i>;-+£ + v,l + 2v ) +2ix^j . (B3) 



wr(| + £ + v) 
Then, the integral I,j can be written as 

Tij = ILeW"-" (v( nx ,u + n 2 , 00) - T>{m,v + n 2 , 0))* (x>(ni, v + h 2 ,oo) - V(n x , v + h 2 , 0)) 

i n 1 +n 2e in(u-u") f fa^fa} (v(n 1 , v + n 2 , 00) - £>(ni, v + n 2 , xi)) , (B4) 



v- 

4 



where (ni, 712), (fii, 712) = (0, 0), (1, — 1), (1, 0) for i,j = 1,2,3, respectively. Here note that T)(£,v,x) gives a finite 
value in the limit x = 00 (see appendix IB II for the derivation) : 



1 ; i^siwv\T{\/2-£-v) T(l/2-£ + v)j y ' 

cos™ 1 J T{£+1) T(\+v)T(\-v) ' 1 j 

where we used an ie- prescription to drop oscillatory terms cx e 2rx at infinity x — > 00. It should be also noted that the 
asymptotic behavior of T>{£, is, x) around x = is given by 

T>(£, u, x) = 2 " T{V) ( x-^- v + 0(xi +£ ~»)) + e -**"-^T^( xh+l+ " + 0(x i+i+ l) , (B7) 
iir(±+£ — v)\ ) iir(j; + £ + v) V / 

which is singular for Re[l/2 + £ — v\ < 0. However, as mentioned in appendix 1X1 this kind of singularities cancel out 
in the calculation of Cy's and we drop them in the following calculation. We therefore rewrite (|B4j) as follows: 



lij = |e^^*)(-^ ni - fil+n2 -" 2 (x'(n 1 ,r/ + ri2,oo))*P(n 1 ,i/ + fi 2 ,cx)) 

POO . 

_ i n 1 +n 2e i v (,- v n I ^J^.( a;i )^2?(n 1 ,i/ + n2,cx>)-X>(ni,i/ + na,xi)J 



1 o-n, -ft, , 1 xn a -n, +n 2 -n 2 (j + ^* + "a)m (f ~ ~ + " + "a)fii (j~ ^ ~ fajgj 



i-iy 



4 cos 2 ™ v ' r(i + m) r(i + ni) 

dxiAi(xi)\V(ni,p + n2,oo) - V{hi,v + h 2 , xi)j , (B8) 



-ni+n 2 \-K(y—v* 

4 



with (a) m = 



F(a + m) 



r(a) • 

We next perform the xi-integral in (|B8[) . To perform this kind of integrals, we use a trick of resummation (22l |. 
Expanding .A, in a; and using the identity for Besscl functions, 

x J v (x)e = ^ a v m x m with a v m = , '' , B9 

' to!v7t1 (m + 2i/+ 1) 

rn=0 v \ / 

Ai is rewritten as 

A = i^ x -h+m e i* . -» (e^^+^J,,^^) - J_ (y , + „ 2) (x)) 

00 00 
™* < r * +,l2 a;-2+™+' l i+™2+^* _ (_!)n 2 ^ a -(^+n a ) a .-|+m+ni-n2-i/*^ j (B10) 

m— m— 
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where again (711,722) = (0,0), (1,-1), (1,0) for i = 1,2,3, respectively. Using this expression, the xi-integral in 
can be written as a sum of integrals in the form 



dx x v (T>(e, u, 00) - V{1, u, x)) . 



(Bll) 



Integrating by parts and using d x D(t, v, x) = A(£, v, x) and x 1+p A(l, v, x) = A(l + I + p, is, x), we rewrite it as 



dxx p (D(l, v, 00) -V(£, v,x)) 



l+p 



{V(t, v,oo)-V(l,v, x)) 



dx—^—A(l+£ + p, u,x) 
l+p 



Jo 



j±- (V(l + * + „, v> 00) + [x 1+ PV(£, u, x) - X»(l + 1 + p, v, x)\ x ^ Q ) 



l+p 



v(i + e+p, 00) 



(B12) 



Here we again dropped the contribution from T> at x = 0, which vanishes or cancels out in our calculation as mentioned 
earlier. Then, the Si-integral can be written as follows: 

_ Z^ni+"2 e f ir(u— v*) 



■l+n 1 +h 1 +n 2 



4 sin ttv* 



dxiAi(x\)[ V(hi, v + n 2 , 00) — T>(h±, v + fi2,xi) 



Ex 



2=0 2 



+ 771 + Til + 77 2 + V 



-V(l/2 + m + Tii + ?i2 + v* + h\,v + fi2, 00) 



e 2 



1 



-(i/*+Ti 2 ) 



2+ 111 + n l- n 2 



0(l/2 + 777 + Ul - 72,2 - V* +77i,^ + 772,Oo)) 

V* / 



(B13) 



After some lengthy calculations, we obtain 



7T I 



l+ni+ni+n 2 



a 



±(i/*+n 2 ) 



4 sin 7W* 



2+771 + 77i ±7l2±fc / 
2~ Til -ni — 2 ^ ^ym+fi2 



-2?(l/2 + 777 + 77l ± 772 ± ^* + 77-1, 7/ + 77 2 , 00 ) 



ShlTW* (I + TO + m ± (l/* + 77 2 )) (| + 777 + (l/* +n 2 )) 1+ni+fti 

f (777 + l)„ 1 ± rl9+fllT f l , (-to - ni + n 2 - 77i + 77 2 + 2i/*)„ lTn , + f ll ± fl9 for real 1/ , 
x < 

[ (to + l)„ 1±Il2+fll±fl2 (-to - 77i + 77 2 - 77i ± 77 2 + 2^*) IllT „ 2+ f ll = F f l2 for pure-imaginary 2/, 



(B14) 



where T(z)r(l - z) 



was used. Finally, it is necessary to re-sum (|B14[) with respect to to. In the case of 



sm7rz 

Re[2n] (rij = hi = 0), for example, we perform the resummation as follows: 



l e i<v-»*) J dxiA 1 (xi)(v(Q,v > oo)-V(Q,v,xifj 
(-l) m e i7TV (-l) m e- im 



4 sin ■kv 

i 

4 sin ~kv 



E 



(m + i + v) 2 (to + | - i/) 5 



[e^$(-l, 2, | + 1/) - e- i7ri/ $(-l, 2, i - i/) 



16 sin7ri/ 



7^ 



2 

3 i/ 
4+2 



■,(D 



1 7/ 

4+2 



16 siwKV 



^"(j-D-^a-D 



i ^ 



Here $(z,s,a) and ip^ n '{z) are the Lerch transcendent and the polygamma function, respectively: 



$(z, s,a) = Y . Z " 1 . , ^W(z) = (-1)" +1 77! V 7 1 -— T 



m=0 



m=0 



(to + z) r 



(B15) 
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which satisfy the following relation: 

#(-l,n + l,z) = (-l)"n!2""- 1 

We then obtain 



(n) 



Z+l 



-V 



(<<) 



(B16) 



Re[Jn] = 



4 COS 2 TTI/ 



Re 



z e 



16 shi7w 



♦ m (! + 5)-* m (l + 5) 



1 V 



i e 



16shi7w 



1 i/> 



(B17) 



which reproduces the result in [22J. In a similar way, we can obtain analytic expression for Xjj's, and the results are 
summarized as follows: 



1 



13 



t-31 



7T>-i|>-!l 2 1 



(e i7r "(2i/ - 3)(2i/ - 5) - e~ ilUJ (2v + l)(2i/ - 1)) 

+ 7^7r4 (2t,-l) 2 (2i/-3) 2 fe^$(-l,2 i + J/)-e-^$(-l,2,f -^)) for real z/ 

256 sin 7W V z z / 



16 cos 2 7Tf 



16 128 sin ?ry 



16 COS 2 7TI/ 
2 



16 128 sin Try 



-3 - 82i/ + 12i/ 2 + 72i/ 3 
3 + 2y 



+ 56j/ 3 ^ 



1 + 2u 



1 sin tti/ 



I33 — 



256 sin7w 



i/(-3 + ^(V™ #(-1, l.f + 1/)- e^^-l, 1, -i - 1/)) 

(4^ 2 -l)(4 I / 2 -9)(e" ,y $(-l,2,| + I /)-e~ i7rl/ $(-l,2,i-^)) for pure-imaginary v, 



16 COS 2 nV 
i 



1 + (e-(2, - 1)(2„ - 3) - e--(2, + 1)(2* + 3)) 



256 shi7w 
7r 2 (^-i)(^-|) 



16 128 sin nv 
{2u + l) 2 (2is - l) 2 (e i7n/ #(-l, 2, § + v) 



e -*™$(-l,2,§-i/)), 



(B18) 



(B19) 



S COS^ nV 



16 SLTLITV 



e mu {2v- I) + e~™ v {2v + 1) 
(-3 + 4*/ 2 ) (e l ™$(-l, \,\ + v)- e-^^i-l, 1, ± - 1/)) 
(-3 + 2i/)(-l + 2i/)(e*"'$(-l, 2, i + 1/) - e-^$(-l, 2, \ - u)) 

(e <w '$(-l 1 l,-j+0 



16 sin7ri/ 



TrV - !)("* - §) , i(-l + 4^ 2 ) 



(B20) 



I COS^ 7TZ/ 



32 siii7Ti/* 



*(-l,l,4-i/ 



i(4i/ 2 - 5) 



f e""$(-l, 1, § + i>*)- e-^^C-l, 1, § - 1/*) 



32 sin tti/ 
i(-3 + 2z/*)(-l + 2^*) 
32 sinm/* 

7^(1/ +!)(«/ -|) 



$(-1,2 i+i/*)- e - 



#(-1,2,1-1/*) 



(B21) 



32 sin7r^ 



+ U nu {2i/ - 1) + e- i7ri/ (2y + 1)) 

16sm7r^V / 

(V^^-l, 1, i + 1/) - e- i7rl/ #(-l, 1,1- i/jj 
■ (>"#(-!, 2, 1 + 1/) - e-^#(-l, 2, 1 - 1/; 



8cos 2 7rz/ ' 16sin7r^ 
t(l + 4i/ 2 ) 

16 sill7Ti/ 

+ 4t/ 2 



(B22) 



+V + \)( V * 1) | z(-l + 4^ 2 ) 
8 cos 2 tti/ 32shi7ri/ 
i(3 + 4z/ 2 ) 



^$(-1, 1 1 + v ) - e - inv <$>(-l, 1, 1 - v) 



32 sinTrz/ 
z(-l + 4i/ 2 ) 
32 sva.ni/ 



(e ? ™$(-l, l.f + i/)- e^^-l, 1,1-1/] 



e Z7ri/ #(-l, 2,| + i/)- e" i7r,/ #(-l, 2, | - v\ 



(B23) 
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T>-il> + |)(f*-f) , i < „ , , .-w-i + 21/ 



23 



f e 4 ™* (-11 + 12^ 2 ) + e- mv * — — (1 - 20^* + 20^)") 

V —3 + 2z/* / 



16cos 2 tw 16 128 sin to/* V -3 + 2^* 

+ 2v*){-l - Av* + Av 2 ) 



-(V™*$(-1, 1, -| + i/*) - e- ?my *$(-l, i/*)) 



32 sin7r^* 

_ })(■,• + !) _ 1 _ + * _ 20i/ 2Qf e _„ _ 24 „ la/ x 

16cos^7w 16 128sin7ri/v 1 + 2v I 



i(-\ + 2v)(-\- Av + Av 2 ) 
32 sin7ri/ 

i (-3 + 2^)(-l + 2^(1 + 2^)/ . 3 



(V™$(-1, 1, | + v) - e-™ v $(-l, 1, § - j/)) 



256 sin 7W 



e 



$(-1, 2, | + i/) - e- m "$(-l, 2,|-i/) . (B25) 
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Here we used (—1)'™ = = — . As we displayed in figure[TJ the analytic results obtained in this appendix 



1 + 1 2 

m=0 

and the numerical results for r, = 1 well coincide with each other 



1. Asymptotic behavior of T>(t,u,x) 

In this subsection we derive the asymptotic behavior (|B5|) in the limit x — > oo of the function "D(£, x): 

2 v x^ +t '~ v T(v) /ll 3 \ 

V{t,vx) = — ^— -i-f 2 F 2 - - i/, - + £ - i/; - + £ - v, 1 - 2v- 2ix) 
i7r(i + £ — v) \2 2 2 / 

1 "- <7n/ 1 V 2 (- + ^-+£ + I ^+£ + M + 2^ a; ) . (B26) 



i7r(i+^ + ^) ' *\2 '2 
We use the following asymptotic expansion of hypergeometric functions 

2 F 2 (a u a 2 -M,b 2 ;z) = If^Jf^ z^-^Y^c^ 

r( ai )r(a 2 ) ^ 



r(6i)r(6 2 )r(a 2 -ax) ,11/^ 
(-z) ^-FiCoijOi - h + l,ai - 6 2 + l;oi - a 2 + l;-l/2) 



r(o 2 )r(6 1 - ai)r(6a - oi) 

r(6 1 )r(6 2 )r(a 1 -a 2 ) 



(—z) a2 3-Fi(a 2 ,a 2 — 61 + 1,02 — 62 + 1; <X2 — ai + 1;— l/z) 



r(ai)r(6i-a 2 )r(62-a 2 ) 

(B27) 

where c^'s are numerical factors independent of z. The first term gives an oscillating term ~ e 2tx , which can be 
dropped by an ie-prescription. Then, let us consider the contribution of the last two terms. They are respectively in 
the form 

x e e -i*(»-m 3Fl (i/2 + u,l/2-v, 1 - £; i/(2x)) 



i\J2Tx£ sin7ri/ 



and 



■ lr J , lml + l 4 r(-Q(-2z)-^-^^- 1 / 2 ) 3 Fi(l/2 + e + v,l/2 + t-v,0;l+t; i/{2x)) , (B28) 
i\/2tt bid. nit 1 (1/2 — t — v) 



— x i e-^ (u - 1/2) 3 F 1 (l/2 + v, 1/2 - v, -I; 1 - i;i/(2x)) 

i\l2-n£ ms.-KV 

7*=7i 1 rn /9 + ^ r(-^)(-20- < e-i^- 1 /2) 3j r l(1/ 2 + £ + „, 1/2 + £ u, 0; 1 + £; i/{2x)) , (B29) 
i\/27r smirv 1(1/2 — £ + 2/J 



28 



where note that hypergeometric functions p F q (ai, . . . , a p ; b\, . . . , b q ; z) are symmetric under the permutations of aj's 
and those of Vs, respectively. The first term in (|B28|) and that in (|B29|) cancel out and we obtain 

V{1, v, x) = , 1 , 1 W£ + \ - iy | r(-^)(^2z)- £ e -^^- 1 /2) 3Fl (i / 2 + I + „, 1/2 + 1 - „, 0; 1 + 1; i/{2x)) 
i\J2n smirv 1 (1/2 — I — v) 



1 1 Y(\j2 + i + v 



Jr(-^)(-2i)-^-^( l '- 1 /2) 3j p 1 ( 1 /2 + 1 + v , 1/2 + 1 - v , 0; 1 + 1; i/(2x)) 



iy/2nsimTPT(l/2-i + v) 

(B30) 

where note that we did not use any approximation so far. Finally, taking the limit x — > oo, we conclude that 
where we used 3^1(01, 02, 03; b; 0) = 1. Using the identity 

r(*)r(i-*) = ^!— , (B32) 

sin ttz 

we can also rewrite (|B31|) as follows: 

0F(l-»)e-i- , 1 T{l + i + v )T{\ + Z-y) 

V(£,u,oo) = (— li) -—- — — —, — —, , (ass) 

v ; costt^ v ; T(£+l) r(| + i/)r(i - v) ' 



which reproduces the result in [22( for I = 0. 
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